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This is a summary of the work accomplished during the second half 
of 1984. Some results from work accomplished in January 1985 are 
i ncluded. 

The main emphasis of the study is to investigate the feasibility 
of rain volume estimation using satellite data following a technique 
recently developed with radar data, called the Area Tire Integral. To 
accomplish this task, case studies were selected (during January 1984 
— see 1st Semi-annual Report) on the basis of existing radar and 
satellite data sets which match in space and time. The radar and 
rapid scan GOES satellite data were collected during the Cooperative 
Convective Precipitation Experiment (CCOPE) and North Dakota Cloud 
Modification Project (NDCMP). 

Four multicell clusters were analyzed in 1984. The small 
mesoscale areas of Austin and Houze (1972) or the mesoscale cloud 
systems described by Super and Heimbach (1980) as cells bigger than a 
cumulus cloud but smaller than a typical mesoscale convective complex 
were considered "clusters." Two clusters were selected on each day, 

12 June and 2 July, as illustrated in Table 1. The 12 June clusters 
occurred during daytime, while the 2 July clusters during nighttime. 

The radar data were on hand at the Institute of Atmospheric 
Sciences, South Dakota School of Mines and Technology. Magnetic tapes 
of digital satellite data were acquired by the Department of 
Atmospheric Sciences, Colorado State University, from the Bureau of 
Reclamation, the University of Wisconsin through the National 
Environmental Satellite Data Information Service (NESDIS), and the 
GOES data archive at CSU. The satellite data were processed at the 
IRIS/DRSES (Interactive Research Imaging System/Direct Readout 
Satellite Earth Station) facility located at CSU. The radar data were 
processed at SDSM&T. Details concerning data processing were 


No. 

Cluster 

Identification 

TABLE 

Date 

1981 

1 

Time Period 
(GMT) 

Satellite 
Data Type 

1 

1A 

12 June 

1627-1925 

infrared, visible 

2 

IB 

12 June 

1716-1925 

infrared, visible 

3 

XI 

2 July 

0035-0453 

inf rared 

4 

X2 

2 July 

0221-0434 

infrared 
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described in the first semi-annual report. Once the satellite data 
were processed into IRIS, the routines for navigation remapping and 
smoothing of satellite images were performed. The visible counts were 
normalized for solar zenith angle. The GOES satellite infrared data 
are calibrated and give count values which can be related to equiva- 
lent black body temperatures. The radar data tapes were processed 
following procedures similar to those described by Schroeder and 
Klazura (1978). Digital printouts of the dBz-values at low-tilt angle 
were prepared and converted from radial to a rectangular coordinate 
system. A radar sector of interest was defined to delineate specific 
radar echo clusters for each radar time throughout the radar echo 
cluster lifetime. The radar sector of interest was used to locate the 
convective cluster responsible for the rainfall in the remapped 
satellite data. The radar echo clusters within the defined sectors of 
interest were identified by drawing a "box' - around the cluster for 
each radar scan. The coordinate of the "boxes" were entered into a 
computer program that calculates the cluster echo areas of >_ 25 dBz 
reflectivity thresholds nnd the corresponding radar ATI. The rain 
volume for each cluster was then computed using an optimized Z-R rela- 
tionship (Smith et al . , 1975). Next, a satellite sector of interest 
was defined by applying small adjustments to the radar sector using a 
manual processing technique. This was done to avoid cloud features 
suspected of not being detected by the radar. It was (and it will be) 
one of the most delicate and difficult tasks of this investigation 
because of the following considerations: 

1) Satellite and radar systems respond to different 
characteristics of clouds at different atmospheric 
levels. During a cluster lifetime, the geometry of the 
cluster as viewed from a low-tilt radar scan and a 
satellite picture of the cloud top changes dramatically. 

2) Time differences between satellite and radar data sets 
can vary by as much as 10 min during which the evolution 
of the convective system continues. 

3) Spacial positioning of the satellite observations 
relative to the radar location necessitate error 
inclusion due to limited accuracy of geometric 
corrections. 

4) Location correction for a cloud as a function of 
height above the earth's surface leads to further 
uncertainties. 

5) Vertical wind shear advects the top of the clouds 
downwind from the location of the radar echo area. 

The direction of vertical wind shear is toward the 
northeast for the two case days. 

6) Another consideration in definition of the matching 
satellite sector of interest involves the inclusion 
or exclusion of cirrus spissatus. If the cirrus 
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debris appeared to be completely detached horizontally 
from the radar echo cluster, the cirrus cloud cover 
was excluded from the satellite sector. A multicell 
cluster may experience reg r, owth during its lifetime. 

This was the case with cell 1A. On the other hand, if 
the radar echo cluster wa r , closely associated with the 
cirrus spissatus, cirrus clouds were included in the 
satellite sector. Cell 1A, 12 June 1981, exhibited 
continual redevelopment on the southern flank. 

The cirrus debris makes definition of the end of a cluster (as a 
convective entity) from the satellite images difficult. The best way 
to avoid such difficulty is to define the maximum development of the 
cluster and consider it as the ending time of the analysis. The ATI 
technique is also valid if only the growing period of the cluster's 
lifetime is considered (Doneaud et a 1 , , 1984). The time of the maxi- 
mum development of the cluster should be determined using both 
satellite and radar data. With further research, it is hoped that 
only satellite data may be used for this purpose. 

The satellite brightness and temperature pixel counts above 
defined thresholds within the satellite sector of interest were counted 
by computer and converted to area (one pixel = 2.47 km 2 ). Histograms 
of area versus brightness counts were produced for each satellite time 
step for each data type for the lifetime of the storm. Graphs are 
then constructed with the ordinate representing digital counts and the 
abscissa being area multipled by time intervals (km 2 x hrs). The 
abscissa represents the satellite ATI equivalent for the whole life- 
time of the storm. A digital count area matching the radar area 
multiplied with the corresponding time interval is also determined. 

Such histograms were produced for the four analyzed clusters 
(Fig. 1-6). 

Different radar and satellite quantities for every time step e or 
the four clusters are listed in Tables 2-5. The step-by-step evolu- 
tion of the echo area multiplied by the time interval and of the 
threshold count value reflect the multicell characteristics of the 
clusters, yet still exhibit similar trends, although there are time 
lags between radar and satellite data. 

Trends of radar and satellite products (from Tables 2-5) are 
displayed in Figs. 7-14. The similarity of the evolution of these 
curves, as well as differences and lags between radar and satellite 
products, are evident. 

Figures 7 and 8 display temporal evolution of satellite and radar 
characteristic quantities for cluster 1A. Two growing and two 
decaying periods are visible in the trends of the radar data empha- 
sizing the multicell character of the cluster. A hint of the first 
growing period is present in the satellite data trend primarily 
because the clusters initial stage was overlaid by debris from pre- 
vious activity. The temporal variation of the radar and satellite 
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Fig. 2: Graph of ATI values for each satellite visible digital count for cluster 1A 

The dashed line represents the digital count for the satellite determined ATI. 






ATI CELL IB INFRARED 1711-1929 GMT 



Fig. 3: Graph of ATI values for each satellite Ik digital counts for cluste" IB. The dashed line 

represents the digital count for the satellite determined ATI. 
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ATI CELL X2 INFRARED 0218-0502 GMT 



Fig. 6: Graph of ATI values for each satellite IR digital counts for cluster X2. The dashed line 

represents the digital count for the satellite determined ATI. 
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13 


increment for cluster 1A. 






Fig. 8: Temporal evolution of maximum echo height, of maximum reflectivity factor (dBz) and of 
satellite count value (°C), matching the interpolated radar area multiplied with time increment 
for cluster 1A. 



characteristic variables for cluster IB is presented in Figs. 9 and 
10. The radar echo area, the IR maximum counts, and the IR counts 
matching radar echo areas evolve similarly (Fig. 9), except for the 
decaying phase of the cluster where the cirrus debris keeps the IR 
counts high. The maximum echo height and the maximum reflectivity 
trends (Fig. 10) are opposite to the IR count matching radar area. 
Radar trends imply that a weak secondary peak is co-located with the 
peak in the IR count trend. Again these differences have to be 
attributed to the multicell nature of cluster IB. The temporal 
evolution of the radar and satellite data for cluster XI are dis- 
played in Figs. 11 and 12. The trends are similar except for the 
initial stage of the cluster. Remnants of previous activity preclude 
initial observations. The cluster X2 (Figs. 13 and 14) demonstrates 
a growing and a decaying period. Its multicell character is self 
evident as the growth and decty of individual cells are indicated. 

Time lags of maximum development between radar and satellite 
variables is also evident. 

The principal goal of this investigation is to compute convective 
rain volumes over stationary or floating target areas by considering 
independent satellite data. As such, the key element is to determine 
ATI from satellite data without using radar returns (as was done 
here). A satellite quantity defining the most appropriate digital 
count matching the radar ATI for a given cluster needs to be found. 

For instance, curves correlating satellite digital count thresholds 
(matching radar ATI's) and satellite IR maximum count values (°C) or 
satellite IR maximum count values (°C) averaged over the duration of 
the cell (or only its growing period) have to be identified. Such 
quantities for the four clusters analyzed are listed in Table 6. It 
is evident that the digital counts matching the radar ATI's are a 
function of cluster characteristics. Similarities between evolution 
of satellite counts matching radar ATI values and satellite IR maximum 
count values do exist; though cell 1A (which was not entirely recorded 
in the radar screen) seems to be an outlier. Certainly only four 
clusters are far from being sufficient for a correlative analysis. 

Work is continuing. Some other 5-6 clusters are scheduled to be 
processed and analyzed in 1985. 


********** 


An extended abstract of the project's progress has been published 
in "Global Scal-t Atmospheric Processes Research Program Review" held 
at the NASA Gocaord Space Flight Center, Greenbelt, Maryland, 

August 8-10, 1984, pp. 203-209. 

A manuscript entitled "Convective Rain Rates and Their Evolution 
During Storms in a Semi -Arid Climate" coauthored by A. A. Doneaud, 

S. Ionescu-Niscov, and J. R. Miller, Jr., has been published in the 
August issue of the Mon . Wea . Rev. , 112 , 8^, 1602-1612. 
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CELL # IB 12 JUNE 1981 
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Increment for cluster IB. 




time increment for cluster IB. 
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Fig. 11: Temporal evolution of radar echo area (km 2 ), of maximum satellite count value (°C) 
and of satellite count value (°C), matching the interpolated radar area multipled with time 
increment for cluster XI. 
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time increment for cluster XI. 


CFLL # X2 2 JULY 1981 
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increment for cluster X2. 
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TABLE 6 
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A paper entitled "An Attempt to Extend the ATI Technique to 
Estimate Convective Rain Volumes Using Satellite Data" coauthored by 
A. A. Doneaud, J. R. Miller, Jr,, L. R. Johnson (SDSM&T), and 
T. H, Vonder Haar and P. Laybe (CSU) was presented in a poster 
session and Included in the preprint volume of the 22nd Conference 
on Radar Meteorology, 10-13 September 1984, Zurich, Switzerland 
(Attachment 1), 

Statements are included in the last papers acknowledging NASA's 
support of the research grant. A list of consulted bibliographic 
references is included in the manuscripts published. 
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